Strange nonchaotic attractors in driven delay-dynamics 

Awadhesh Prasad^ , Manish Agrawalt , and Ramakrishna Ramaswamy * 



^ Department of Physics and Astrophysics, University of Delhi, Delhi 110 007, India 
i School of Physical Sciences, Jawaharlal Nehru University, New Delhi 1 10 067, India 

Abstract: 

Strange nonchaotic attractors (SNAs) are observed in quasiperiodically driven 
time-delay systems. Since the largest Lyapunov exponent is nonpositive, trajectories 
in two such identical but distinct systems show the property of p/ioie-synchronization. 
Our results are illustrated in the model SQUID and Rossler oscillator systems. 



1 Introduction 

Extensive studies over the past twenty or so years on strange nonchaotic attractors |[T] 
have helped to establish that such behavior is not exceptional dOIU. By now several 
examples are known where it may be clearly shown that the dynamics is both strange 
and nonchaotic LL^^ TJ, namely that the attractor has a fractal geometry, and that the 
largest Lyapunov exponent is either zero or negative. A number of questions relating 
to the origins of such dynamics remain open. Most studies (analytical or numerical) 
of SNA dynamics have been on low-dimensional quasiperiodically driven systems. 
One open question relates to the need for external forcing. All known examples of 
systems where SNAs occur have a skew-product structure, and it is still not entirely 
clear whether this feature is necessary or merely sufficient. Further, in all known 
examples the external drive is quasiperiodic in time: it is also not clear whether this 
form of drive is necessary for the creation of strange nonchaotic motion or merely 
sufficient [3|. 

For chaotic attractors, there can in principle be several positive Lyapunov expo- 
nents (the phenomenon of hyperchaos). Although the SNAs are geometrically similar 
to chaotic attractors, any further analogy between them is clearly not possible, and 
thus it is of interest to investigate the nature of SNAs in high-dimensional dynamical 
systems. 

In the present work we consider model time-delay quasiperiodically driven dy- 
namical systems and study the nature of the dynamics. The inclusion of time-delay 
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through a diffusive self-feedback coupling term makes the system effectively infinite- 
dimensional. The SNAs that are created, however, continue to be low-dimensional 
IISl and are very similar to those created in the absence of time-delay seen in earlier 
studies f9l. A second motivation for the present study is practical. In order to create 
SNAs in an experiment, the physical set-up may involve time-delays and thus it is 
important to know the effects of including such coupling. Time-delay coupling has 
attracted considerable interest since this can cause interesting dynamical phenomena 
such as amplitude death UOJ or novel bifurcations lfT2l[TlJ . 

Since all Lyapunov exponents on SNAs are nonpositive, a distinctive property 
of such attractors is that trajectories with different initial conditions do not separate 
from each other. Indeed, on two identical but separate systems, trajectories started 
with the same phase will completely synchronize with each other ITSl . If the ini- 
tial phases are different, then trajectories show phase-synchronization 11411 . Similarly, 
on a given SNA, trajectories starting from different initial conditions also coincide or 
phase-synchronize in this manner, which is one of our aim in this paper 

In the next section of this paper we study the occurrence of SNAs in time-delay 
driven dynamical systems. This is followed by a discussion of phase synchronization 
in Section III. Finally, a brief summary is given in Section IV. 

2 Delay SNAs 

We study two model quasiperiodically modulated systems. Among the earliest demon- 
strations of dynamical systems with SNA was the driven and damped pendulum equa- 
tion 19] used to model a driven SQUID with inertia and damping. We introduce a delay 
feedback term, to get the equation of motion as 

X + "fx + X + q2Sm2'KX = gi sin(a-'it) + /?sin(w2t) 

+e[x{t - t) - x{t)]. (1) 

The drive is made quasiperiodic by requiring the ratio of the frequencies to be an 
irrational number Here we take this ratio as inverse golden mean, i.e. uji/lu-? = uj = 
(75-1)72. 

In a similar manner, we also consider a driven Rossler system flSl with delay 
feedback and quasiperiodic parameter modulation, 

X = — y — z 

y — X + a{l + -{cost + cos LLjt))y 

+e[y{t-T)-y{t)] 
z = 0.1 + z(.x- 14). (2) 

The above dynamical systems are integrated using standard techniques |fT6ll and we 
compute the largest several Lyapunov exponents 1 17| as a function of parameters. For 
a wide range of parameter values, we find that the dynamics is on nonchaotic attractors 
in both systems, and representative results are shown in Fig. 1. 
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The largest Lyapunov exponent is shown in Figs. 1 (a) and (b) for the SQUID 
and Rossler systems respectively. (See the caption for details of the parameter values). 
The symbols T and C denote quasiperiodic torus and chaotic attractors respectively, 
and AD indicates the amplitude-death region, where oscillations die to a fixed point 
lfT2l . The corresponding variance of finite-time Lyapunov exponents, which has been 
used as an order parameter for detecting the torus to SNA transition (see Ref. ifTSl 
for details), is shown in Figs. 1 (c) and (d) and this suggests the presence of SNAs 
in the region marked by arrows. Examples of such nonchaotic attractors are shown in 
Figs. 1 (g, h) for the two cases. The route to SNAs appears to be via fractaUzation 
EHU (see Figs. 1 (e, f)); in this transition, a torus attractor gets wrinkled progressively 
as parameters are increased and transform to SNAs y_9| with a slow increase in the 
variance of finite-time Lyapunov exponents as shown in Fig. 1 (c, d). 

A number of measures 1201 [3] |2T] have been suggested in order to quantitatively 
confirm that the attractors are indeed SNAs. We compute the partial Fourier sum ||2TI . 



where fl is proportional to the irrational driving frequency uj fTP\ and is the time 
series of one of the dynamical variables. The graph of Re T vs. Im T gives a "walk" on 
the plane, with mean square displacement|r(ri, N)\'^. The singular-continuous nature 
of the SNA specti-um imphes |20| the scaHng \T{n, 7V)p ^ iV^ with 1 < ^ < 2. Plots 
of log |T(ri, iV)p versus log N in Figs. 1 (i,j) for the SNA dynamics corresponding to 
Figs. 1 (g,h), show linear behavior, with slopes fi « 1.67 (for the SQUID system) and 
fi « 1.59 (Rossler system), suggesting that both the attractors are indeed strange and 
the dynamics is nonchaotic. 

3 Synchronization 

Since the largest Lyapunov exponent is negative, trajectories on SNAs have the prop- 
erty that they eventually coincide and become identical (typically with a time-shift). 
Explicitly, for two different initial conditions denoted x(0) and x'(0), the distance 



for some t. When it is possible to define a phase variable, then this is similar to a 
phase-shift. 

Consider now two identical systems where the dynamics is on SNAs. Trajectories 
starting from arbitrary initial conditions will therefore show synchronization. If the 
initial phases are the same, the synchronization will be complete ITSll . but more gen- 
erally, this will be a phase-synchronization 1 14] . Since this phenomenon occurs in the 
absence of any coupling between the two systems, this notion of synchronization also 
applies to two trajectories in the same system. 

The distance between two trajectories with random initial conditions is shown as 
a function of time in Fig. 2 (a) (the solid line) for the SQUID system. This quantity 
rapidly goes to zero. For the Rossler case, earlier work has shown that it is possible 
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X(t)-X'(t+T)|| ^0 
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Figure 1 : Results for the SQUID system with parameters r = 2, coupling strength 
e = 0.01, qi = 2.768, g2 = 7 = 2 and the Rossler system with t = 0.9. The 
Lyapunov exponents are plotted as a function of drive amplitude f3 for the SQUID 
system in (a) and as a function of a in the Rossler system in (b). The variance in the 
Lyapunov exponents are shown in (c) & (d) for the cases (a) & (b) respectively. The 
representative Poincare sections of periodic torus and corresponding SNAs are shown 
for the two systems in (e) /3 = 1 (y = i), & (f) a = 0.1, and (g) /3 = 1.6 & (h) 
a = 0.114. Here the Poincare section are taken at x = 0, and then 9 = mod{t, 2tt) 
are considered. Shown in (i) and (k) are the "walk" displacements, T{il, N) vs N on 
a logarithmic scale for SNAs corresponding to (g) and (h) respectively at = uj / 4. 

to define a "phase" as (j) ^ ta,n~^ {y / x) fl4]. For two trajectories on SNAs, the phase 
difference A(f> = 4>i — (j)2 = Awi does not grow with time (see Fig. 2 (b)) unlike 
what happens in the case of very morphologically similar chaotic Rossler attractors. 
As shown in Figs. 2 (c)-(d) on SNAs (namely in the parameter range indicated by 
arrows ) the motion is phase locked in both the SQUID and Rossler systems. 

4 Summary 

In the present work we have studied representative nonlinear time-delay dynamical 
systems with quasiperiodic forcing and observe that the dynamics can be on strange 
nonchaotic attractors. These attractors are low-dimensional and trajectories on these 
SNAs have the property of phase-synchronization. We have verified that these features 
are shared by other delay systems such as the Mackey-Glass equations Il22ll . or the 
time-delay Duffing oscillator when subject to an external quasiperiodic drive. 
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Figure 2: (a) Distance between two trajectories starting with different initial condi- 
tions for SQUID, Eq. ^ (solid line) and Rossler, Eq. ^ (dashed line) in SNA region 
at parameters /? = 1.6 and a = 0.114 respectively, (b) The phase difference for 
a chaotic attractor (solid line) at a = 0.2 and SNA (dashed line) at a = 0.114 for 
Rossler oscillator, (c) and (d) show the frequency differences for SQUID and Rossler 
systems respectively. Arrows indicates the regimes of SNAs. The curves in (c)-(d) are 
averaged over 50 random initial conditions sampled uniformly. 
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